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Abstract 

Real-time neutrino oscillation experiments such as Super-Kamiokande (SK), the Sudbury Neu- 
trino Observatory (SNO), the Kamioka Liquid scintillator Anti-Neutrino Detector (KamLAND), 
and Borexino, can detect time variations of the neutrino signal, provided that the statistics is 
sufficiently high. We quantify this statement by means of a simple unbinned test, whose sensitivity 
depends on the variance of the signal in the time domain, as well as on the total number of signal 
and background events. The test allows a unified discussion of the statistical uncertainties affecting 
current or future measurements of eccentricity-induced variations and of day-night asymmetries 
(in SK, SNO, and Borexino), as well as of reactor power variations (in KamLAND). 

PACS numbers: 14.60.Pq, 13.15.+g, 26.65.+t, 95.55.Vj 
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I. INTRODUCTION 



Among the experiments that have estabhshed the phenomenon of neutrino flavor os- 
cillations (see, e.g., [l| for a recent review), several ones are able to accurately tag each 
neutrino event in the time domain. In particular, the Super-Kamiokande (SK) experiment 
[2] and the Sudbury Neutrino Observatory (SNO) can detect solar neutrino events in 
real time through the Cherenkov technique, while the Kamioka Liquid scintillator Anti- 
Neutrino Detector (KamLAND) Q can detect reactor neutrino events in real time through 
the scintillation light. The Borexino experiment (in construction) 5] will also allow real-time 
observations of solar neutrino events through a liquid scintillator technique. 

In all such experiments, the separation of the global neutrino data sample into signal 
(S) and background (B) cannot be performed on an event-by-event basis, but only on a 
statistical basis, providing the total number of events in each class: 

N = Ns + Nb. (1) 
Therefore, within a generic time sequence of events, 

=1,...,Ns+Nb , (2) 

the i-th one can be either a signal event or background event, with probability Ns/N or 
Nb/N, respectively. Testing time variations of the signal is then made difficult by unavoid- 
able statistical fluctuations of the background around its average level. This problem is 
often tackled by binning the events in time intervals At, which should be shorter than the 
signal variation timescale, but also long enough to allow a reasonable statistical subtraction 
of the background. However, as it was emphasized in 0, Q , one can also directly use the 
experimental information in Eqs. and without binning.^ 

In this paper, building upon Refs. 2} , we discuss the statistical significance of possible 
time variations of signals in real-time experiments, by means of a simple unbinned test based 
on the time sequence of events {ti}. In Sec. II we introduce the test. In Sec. Ill we apply 
it to the discussion of seasonal variations induced by the Earth's orbital eccentricity on 
the solar neutrino signal in SK, SNO, and Borexino. In Sec. Ill we discuss the statistical 
significance of possible day-night variations induced by Earth matter effects on the solar 
neutrino signal in SK and SNO. In Sec. IV we discuss the possibility of tracking the time 
variation of the total reactor neutrino flux in KamLAND, with and without the background 
from geoneutrinos or from a speculative georeactor. Finally, we draw our conclusions in 
Sec. V. Technical details about the statistical power of the test and about its comparison 
with the Kolmogorov-Smirnov unbinned test are given in Appendix A and B, respectively. 



II. AN UNBINNED TEST OF TIME VARIATIONS 

Let us consider a real-time neutrino experiment characterized by an event rate 

Rit) = Sit) +B , (3) 



The works |6|, l7| were focussed on the Fourier analysis of time variations induced by vacuum oscillation 
solutions to the solar neutrino problem, which are currently ruled out. 
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where S(t) is a generic time-dependent signal and B is the background, assumed to be 
constant (up to statistical fluctuations). If the function S{t) is known (e.g., from theoretical 
predictions), it can always be cast in the form: 

S{t) = s(l + ^f{t)^ , (4) 

where 5* and (j| are the average value and the variance of S{t) over the detection time 
interval T: 

S = ^j^dtS{t), (5) 

4 = ^[ dt[S{t)-S]\ (6) 

so that f{t) has zero mean value and unit variance in the same interval [0, T]. It is natural 
to ask then how well can one prove experimentally that the signal varies in time, i.e., that 
S{t) is different from its average S or, equivalently, that f{t) ^ 0. 

A simple "unbinned" test of the hypothesis /(t) 7^ consists in showing that the quantity 

Ns+Nb 

/?= E /(^O, (7) 

i=l 

defined in terms of the real-time event sequence in Eq. (j21), is statistically different from zero. 
If the variance (t| can be estimated, then the difference of (3 from zero can be expressed in 
terms of standard deviations ngj, 

Usd = ■ (8) 

We show in Appendix A that, for as/S -C 1, it is 

f3±apc:^Ns^±^Ns + NB , (9) 
which implies the following estimate for ngd'- 

Equation (|l(Jj) provides a useful test of the sensitivity to time variations. It express the 
statistical sensitivity (number of standard deviations rigd) in terms of a theoretical quantity 
(the ratio (Ts/S for the expected signal) and of two experimental quantities (the total number 
of signal and background events, Ns and A''^), independently of specific detector details.^ 
The proportionality to as/S is clearly expected, but the dependence on Ns / a/ Ns + Nb is not 
entirely trivial. A similar dependence arises in the determination of the Fourier components 
of a periodic neutrino signal . In addition, we show in Appendix B that such dependence 



^ Detector details are relevant for the evaluation of experimental systematic errors, including systematic 
variations of the background and of the detection lifetime or efficiency (which can be relevant in some 
cases not considered here, e.g., short-period time variations Q)- Since such evaluation is beyond our 
expertise, we refer to statistical errors only in this paper. 
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arises also in the Kolmogorov-Smirnov unbinned test, and can thus be considered as a rather 
general scaling rule for the statistical sensitivity to time variations. 

Despite its simplicity, Eq. (fTUj) can be useful to understand several issues related to 
the experimental detection of time-varying neutrino signals, without performing numerical 
simulations or data fitting. Representative cases will be discussed in the next Sections. The 
reader is referred to Appendix A for the derivation of Eq. (jU)), upon which Eq. (jlUj) is based. 



III. TESTING ECCENTRICITY EFFECTS 

Within the current large-mixing-angle (LMA) solution to the solar neutrino problem P| , 
the neutrino oscillation phase at the Earth is completely averaged out (see, e.g., and the 
only effect of the Earth's orbit eccentricity {e = 0.0167) is a variation of the signal with the 
inverse square distance, 

S{t) = S (^l + 2ecos'^^ +0{e^) , (11) 

where T = 1 year and t = at perihelion. In this case it is cts/S = V2e, so that Eq. (fTUj) 
can be written as 

n,d = V2£\ — . (12) 

The above equation will be used to test the statistical sensitivity to eccentricity effects in 
SK, SNO, and Borexino. 



A. Application to SK 

In its first phase of operation (SK-I), the Super-Kamiokande experiment has collected 

Ns ~ 22400 (13) 



signal events induced by solar neutrinos P, lld|. The corresponding signal-to-background 
ratio can be estimated through the well-known plot showing the distribution of all events 
as a function of cos6'sun In this plot the background is flat in cos6'sun5 while the signal 
is peaked at cos6'sun = 1 and extends down to cos^gun — 0.5. By adopting the reasonable 
cut cos6'sun ~ 0.5, we graphically estimate from [l3| that 

Nb/Ns - 2.5 . (14) 

We thus expect From Eq. (fT^ a sensitivity to eccentricity effects equal to (in unit of standard 
deviation): 

nsd ^ 1.9 . (15) 

Our estimate Usd — 1.9 might seem too low as compared with the SK-I official analysis 
jlol |. where evidence for eccentricity effects is claimed at ~3(T level. However, we note that 
the best-fit eccentricity measured by SK is a factor ~ 1.5 greater than the true one [lO |. 
This accidental upward shift (^sk — 1-5 e) leads to a corresponding sensitivity enhancement, 
— 1-9 X 1.5 = 2.8, not far from the ~3(T detection level in Ref. 3|. Therefore, it may 



be argued that the intrinsic sensitivity of the SK-I data sample to eccentricity effects is at 
level of ~ 2(7, and that the ~ 3cr evidence in Ref. may partly be due to an accidental 
upward fluctuation of the signal. 
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B. Application to SNO 



The Sudbury Neutrino Observatory experiment has reported physics results for two main 



phases: SNO I (pure D2O, 306.4 hve days) jjj, iJJ, [JJ and SNO II (D2O plus salt, 254.2 



live days) [15|- The total number of solar neutrino events for SNO I+II is 

Ns ~ 5662 , (16) 

including all charged current (CC), neutral current (NC), and elastic scattering (ES) events. 
The background is very efficiently rejected {Nb — 208 events in phases I+II), leading to 

Nb/Ns - 3.7 X 10-^ . (17) 

Our estimate for the SNO I+II sensitivity to eccentricity effects is then [Eq. (jl2j) ]: 

n,d ^ 1.74 . (18) 

Therefore, we argue that the SNO experiment should already be able to see first indications 
for eccentricity effects at > 90% CL.^ If the Nb/Ns ratio is kept at the low current value, 
a future increase of the statistics by a factor of two (three) should bring the sensitivity to 
the interesting level of about 2.5a {3a) in the future. 

C. Application to Borexino 

The Borexino solar neutrino experiment (in construction) is expected to detect, in the 
absence of oscillations, about 55 events per day in the energy window [0.25,0.80] MeV @. 
In the presence of oscillations, the rate suppression depends on the values of the dominant 
oscillation parameters (5m^, 612). Figure 1 shows iso lines of the the suppression factor in the 
mass-mixin g pl ane, superposed to the current LMA solution to the solar neutrino problem (as 
taken from jl6[)- From this figure we derive that the LMA-oscillated event rate in Borexino 
should be roughly equal to 13000 events per year, give or take one thousand events. After 
three years of data taking, the number of signal events should then be around 

Ns ~ 3.9 X 10^ (3 years) . (19) 

One of the goals of the Borexino experiment is to reach a background event rate compa- 
rable or smaller than the signal rate j5|, 

Nb/Ns < 1 . (20) 

If such goal is achieved, from Eq. (fT^ we estimate that, after three years of data taking, 
the sensitivity to eccentricity effects should be definitely larger than 3a: 

> 3.3 . (21) 

Such "eccentricity test" will convincingly prove that the Borexino signal, despite the lack 
of directionality and the nonnegligible background-to-signal ratio, does come from solar 
neutrinos. 



For one degree of freedom, 90% C.L. corresponds to 1.64(7. 
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IV. TESTING DAY-NIGHT ASYMMETRIES 



Within the current LMA solution to the solar neutrino problem (see Fig. 1), Earth matter 
effects are expected to induce r aarg inal day-night variations of the signal S in the SK 
and SNO experiments (see, e.g., [ol. llOl JlS^). Such effects, which decrease with increasing 
Sw? / E^, (where E^, is the neutrino energy), should be null in Borexino. 

Earth matter effects are usually parametrized in terms of a day-night asymmetry A, 

^{Sn + Sd) 

where Sn and Sd are the average signal event rates during night and day, respectively, with 
S ={Sd + Sn)/2. Within the LMA solution it is A > 0. 

For our purposes, we can assume an approximate step-like variation of the signal^ as a 
function of the zenith angle 9^ (a variable more useful than the time t in this context): 

^(nr^^n ^ ^ / -^(^ ~ , cos^^ G [-1,0] (daytime) , . 

^^cost/,j - \ 5(1 + A/2) , COS0, G [0,+l] (nighttime) 

Within such approximation, the signal variance is simply given by 

Using Eq. ()10|). the statistical significance of a day-night asymmetry A is, in units of standard 
deviations, 

A I Ws , , 

= 2il + Ns/Ns ■ (''^ 

We shall use the above equation in a slightly different version, in order to estimate cr^, 
namely, the la statistical uncertainty of A. This uncertainty is obtained by setting n^d = 1 
and A = (Ta in the above equation. The result. 



'I + JVb/JVs 



will be used for numerical estimates in SK, SNO, and Borexino. 



A. Application to SK 

As in Sec. II A, we take for SK-I the values Ns ~ 22, 400 and Nb/Ns ~ 2.5. Our estimate 
for the la uncertainty of the day-night asymmetry A is then 

aA ^ 2.5 X 10"^ , (27) 



^ Figure 2 in 10] shows, e.g., that the step-hke approximation in SK fails only for cos9z ^ 0.4, which 
corresponds to a relatively small fraction (^^1/5) of the total signal. 
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TABLE I: SNO data from phase I (306.4 live days): Total number of events for each of the three 
classes CC, NC, ES and for the total solar neutrino signal S (2nd column), together with their fe 
and u^r components for (^^^/(pe = 1-94 (3rd and 4th column, respectively). See the text for details. 



Class 


Total events 


component 


f^r component 


X 


Nx 




N^/ 


CC 


1967.7 


1967.7 





NC 


576.5 


196.3 


380.2 


ES 


263.3 


203.0 


60.6 


S 


2807.8 


2367.0 


440.8 



which is in reasonable agreement with the value 2.0 x 10~^ quoted by the SK Collaboration 
in Ref. (iq| for the "old" analysis method, based on the comparison of integrated daytime 
and nighttime samples. 

In the same article [l^, the SK collaboration also reported that, by using an extended 
maximum likelihood method, the statistical uncertainty of the day-night asymmetry can 
be reduced from 2.0 x 10~^ to 1.6 x 10~^, the latter being adopted as official SK value 
for aA- We interpret this reduction as an effective improvement in background rejection 
achieved through the SK maximum likelihood method (i.e., lower Nb/Ns and thus lower 



(Ta oc yl + Nb/Ns). We cannot further elaborate on this interpretation for lack of relevant 
information (the SK maximum likelihood analysis in is currently not reproducible outside 
the Collaboration); public release of such information would be beneficial to improve current 
global analyses of solar neutrino data, and to allow combined SK+SNO tests of day-night 
asymmetry effects (see also Sec. IV C). 

B. Application to SNO (phase I) 

In its first phase of operation [l^, l3|, the SNO detector has collected a total of Ns = 
2807.8 solar neutrino events with small background {Nb/Ns = 4.4 x 10~^, neglected in the 
following). The signal events are statistically separated into three classes (CC, NC, and ES), 
each class containing different contributions of events induced by z/g and u^^r- Table I shows 
the content of each class, as obtained by fixing the neutrino flux ratio 0^,T/0e = 3.41/1.76 = 
1.94 and the cross-section ratio a^/ae = 1/6.48 13l|.^ 

The SNO-I result for the day-night asymmetry of the z/g component of the total signal 

(^e) is H 

A = [7.0 ± 4.9(stat.)+}:^(syst.)] x 10"^ , (28) 

under the standard assumption of no asymmetry {A = 0) for the total neutrino flux = 
0e+0/xr- In this section we reproduce with good accuracy the SNO-I statistical error estimate 
= 4.9 X 10~^. On the basis of this successful check, in the next Section we will try to 
estimate the value of after the SNO I-I-II phases. 



^ Variations of the flux ratio (j>^.T/4'e — 1-94 within experimental uncertainties do not affect appreciably 
our results. 
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We remind that the and z/^,- components of the solar neutrino signal in SNO can 
be separated only on a statistical basis. Therefore, the v^. asymmetry A^. (which is not 
directly observable) must be linked to the day-night variation of the total signal S (which 
is measurable), as described in the following. 

The constraint of no day-night change for the total flux (0e + 0/^t) implies that the 
(approximately step- like) day-night variation of 0e, 

A 



1±y) ' (29) 



is associated to the following day-night variation of 



■'ilT ^ ilT 



(30) 



2 



By applying the above variation factors to the and A^^^ components of the total number 
of solar neutrino events Ns-, one gets the global day-night variation 

Ns^Ns^^(n,-N^,^ , (31) 



2 



which, using the numbers in Table I (last row), provides the desired link between Af^ and 
the total signal variation, 

iV5-^iV5(^l±0.76^) . (32) 
The analogous of Eq. ^ is then (for Ns - 2808 and NbI^s < 1): 



in good agreement with the official SNO-I statistical error [l^, also reported in Eq. (j2HI)- 
On the basis of this successful test, we "predict" in the next section the statistical error a^^ 
expected from the day-night analysis of SNO I+II data. 



C. Application to SNO (phase 



In its second phase of operation llq, the SNO detector has increased the solar neutrino 
statistics by 2854.1 events, including 3307.3 CC, 1920.7 NC, and 433.9 ES events. Adding 
these events to those of phase I, one gets the total numbers reported in Table II (where, for 
the sake of simplicity, we have used the same ^/^r/^e ratio as for the SNO phase I). 

By inserting in Eq. (jHT|) the values of Ng and A^^^^ from Table II, we get the following 
day-night variation for the total solar neutrino signal in SNO I-I-II [Ns = 5661.9): 

Ns^Ns(l± 0.63^) , (34) 



so that 



2 



cja ^ ^ = 4.3 X 10"^ . (35) 

^ 0.63 x/iV; ^ ^ 
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TABLE II: SNO data from phases I+II (306.4 + 254.2 live days): Total number of events for each 
of the three classes CC, NC, ES and for the total solar neutrino signal S (2nd column), together 
with their z^g and v^r components (3rd and 4th column, respectively). See the text for details. 



Class 


Total events 


component 


f^r component 


X 


Nx 




N^/ 


CC 


3307.3 


3307.3 





NC 


1920.7 


653.9 


1266.8 


ES 


433.9 


334.2 


99.7 


S 


5661.9 


4295.4 


1366.5 



Therefore, we estimate that the statistical error of the day-night asymmetry Ae for 
SNO I+II should be about 4.3 x 10~^, i.e., not much smaller than for phase I only [see 
Eq. (I2H1)]- Notice that, if the CC sample were hypothetically isolated on an event-by-event 
basis with no background (an irrealistic goal), the corresponding statistical uncertainty of the 
Ue asymmetry would be ctacc = 2/ a/ Nqc = 3.5 x 10~^ for SNO I+II. These considerations 
suggest that, in any case, the uncertainty of the z/g day-night asymmetry in SNO I+II cannot 
be smaller than ~4 x 10~^. This statement will be checked soon, since the SNO I+II official 
day-night is currently being finalized jl9l|.^ 

Can the combination of the day-night asymmetries in SK (as reported in Ref. and 
in SNO I+II [our estimated error in Eq. p5|l ] improve the current determination of the 
LMA oscillation parameters? The answer can be derived from Fig. 2, where we map the 
current LMA region at 99.73% C.L. (as taken from U^) onto the plane charted by the SK 
and SNO day-night asymmetries.^ There is a stron g po sitive correlations between the SK 



and SNO asymmetries, as it was emphasized in Ref. j20| and subsequently in Ref. [21[ . The 
asymmetry rapidly decreases for increasing values of the neutrino squared mass difference 
5m^, for which three representative values are shown in Fig. 2 (including the LMA best fit, 
= 7.2 X 10~^ eV^ (l^). On top of the LMA region, we superpose the SK day-night 
measurement at ilcr from |lfl|, plus a prospective SNO I+II measurement characterized by 
our estimate for the statistical error [Eq. (j35|) and no systematic error] , and by the "luckiest" 
central value (on top of the LMA best fit point). It can be seen that the SK+SNO "error 
box" at ±l(j is larger than the current LMA region at 99.73% C.L. Therefore, although the 
combination of latest day-night asymmetry datum from SK and of the prospective one 
from SNO I+II (our estimate) can provide a useful consistency check of the LMA predictions, 
we do not expect that such data can significantly reduce the current LMA parameter region. 

^ In the final analysis, the current SNO lifetime for phase II (254.2 days 0|) might include additional 150 
live days i.e., about 60% more statistics. In this case, we estimate that the SNO I+II error aA should 
decrease from 4.3 x 10"^ j^q. ^] to 4.0 x 10"^ 

^ The "substructures" of the LMA region in Fig. 2 represent graphical artifacts, due to mapping of a finite 
number of points. 
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FIG. 2: Correlation plot of the day- night asymmetry in SK and in SNO. The slanted region maps the LMA 
solution in Fig. 1 (at 99.73% C.L.), which strongly correlates the predictions for the day-night asymmetry 
in SK and SNO. Such predictions depend mainly on Sm^ (three representative values being shown in the 
figure). The datum with vertical errors bars (ilcr stat.-|-syst.) corresponds to the SK day-night measurement 
from [l^. The "datum" with horizontal error bars implements our estimate for the day- night asymmetry 
uncertainty (±lcr stat. only) from prospective SNO I-fll data (306.4 -I- 254.2 live days), with hypothetical 
best-fit central value. 



D. Application to Borexino 

In Borexino, one expects no day-night asymmetry {A = 0) within the LMA solution 
to the solar neutrino problem. Using Eq. ()26p and the same input values as in Sec. II C 
{Ns = 3.9 X 10^ and Ns/Nb ^ 1), we estimate the accuracy of this null result to be at least 



A±aA= (0.0 ± 1.4) X 10" 



(36) 



after three years of data taking. A check of the null result {A = 0) with percent accuracy 
will provide a useful test that the Borexino detector works as expected, although it will not 
improve the determination of the neutrino oscillation parameters in the LMA region. 
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Time evolution of Reactor Fission Flux at KamLAND 

7 r — \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ r 



D 




Time (day) 



FIG. 3: Variations of the total reactor fission flux power in KamLAND during 540 days of operation (as 
taken from The horizontal dashed line represents the average flux in the whole period. Published 

KamLAND data ^ refer to the first 216 days. 



V. TESTING REACTOR POWER VARIATIONS IN KAMLAND 

The KamLAND experiment 0,113 is collecting events induced by produced in (mainly) 
Japanese reactors. The reactor power demand (and thus the neutrino flux) in a given 
country generally follows a seasonal trend On top of this trend, a strong tempo rary 
reduction of Japanese reactors' power has occurred at the end of 2002 and in 2003 |24| . 
Figure 3 (graphically reduced from jl^]) shows the total reactor flssion flux during 540 days 
of KamLAND operation (starting on March 4, 2002), where the flrst 216 days correspond 
to the flrst phase with published results [22]. In the following, we investigate whether such 
flux variations are detectable with 540-day statistics,^ with or without background, by using 

Eq. (uni). 



We stop at 540 days since, to our knowledge, the reactor flux at later times has not been publicly presented 
by the KamLAND Collaboration so far. 
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A. Case without background 



In the first 216 days, the KamLAND experiment has collected 54 events above the 2.6 
MeV analysis threshold, which cuts away basically all backgrounds ji^. If we naively rescale 
these events using the time variation curve in Fig. 3, we estimate that the statistics should 
be approximately doubled after 540 days,^ 

A^s ^ 2 X 54 = 108 . (37) 

From the same curve in Fig. 3, we obtain that 

^ ~ 0.22 (38) 

over the whole period. Therefore, we estimate that the KamLAND experiment should see 
an indication for reactor flux variations with a statistical significance equal to 

nsd = y - 2.3 (39) 
standard deviations, after 540 days of operation. 



B. Case with background from geoneutrinos or from a georeactor 

Below the 2.6 analysis threshold, the main irreducible known background in KamLAND 
is due to geoneutrinos jl^l. Since geo-i^e background predictions are uncertain jl^, we adopt 
for reference the KamLAND best-fit estimate of Nb — 9 geoneutrino events, to be compared 
with the total reactor signal with no threshold, Ns — 71 events 2^ . 



The previous numbers refer to the first 216 days. At the end of 540 days of operation, 
the total reactor signal should be about doubled (as noted previously), 

Ns ~ 142 , (40) 
while the constant geoneutrino background should be rescaled by a factor 540/216 = 2.5, 

Nb ^ 22.5 . (41) 
Inserting these numbers [and Eq. in Eq. (jTU)), one gets 

nsd ^ 2.4 , (42) 

which is comparable to the no-background estimate in Eq. (j39|) . Reactor flux variations can 
thus be tested at > 2cr (after 540 days) also in the presence of geo-i^e background. 

While geoneutrinos represent a "guaranteed" z/g background, a natural reactor source 
in the Earth's core ("georeactor") [26| represents a more speculative hypothesis, which is not 
endorsed by standard geochemical Earth models (see, e.g., [1^). For the sake of curiosity, 
we estimate that — as a rule of thumb — a georeactor of power p (in units of TW) should 



Of course, a more detailed estimated should take into account power variations and oscillation effects for 
each reactor. 
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add a fraction of about p% to the average KamLAND signal from man-made reactors. 
Assuming then a maximum power p ~ 10 TW [26^], the additional georeactor background 
{Nb/Ns — 0.1) should degrade the statistical significance estimated in Eqs. (jH^ and (jl^ by 
~0.1(j or less. We conclude that, with or without geoneutrino or georeactor backgrounds, 
the KamLAND experiment should definitely be able to see time variations of the reactor 
signal at a level ^ 2a (after 540 days of operation). This check will provide further confidence 
on the origin of the signal from reactor power plants. 

If time variations of the reactor neutrino signal are successfully tested in KamLAND, 
one could, in principle, turn the test around, and try to constrain the amplitude of constant 
backgrounds such as geoneutrino or georeactor events. This kind of tests, together with 
detailed analyses of energy spectra, might improve the discrimination of man-made reactor 
signals from geo-signals. Such analyses could be performed also outside the Collaboration, 
provided that: (1) the neutrino flux "history" of each single reactor is publicly released, as 
it was emphasized in [i^; (2) the energy and time tag of each event (background and signal) 
is also released. The second condition is particularly important to perform unbinned data 
analyses in the geoneutrino energy window, where the few events being collected contain 
very precious information for Earth sciences. We hope that the KamLAND Collaboration 
will take these two desiderata into account. 

VI. SUMMARY AND CONCLUSIONS 

Real-time neutrino oscillation experiments such as SK, SNO, Borexino, and KamLAND, 
can detect time variations of the incoming neutrino flux. We have introduced and discussed a 
possible test for time variations, which does not require time-binning of the events. The test 
provides a simple estimate for the statistical sensitivity to time variations, in terms of the 
signal variance and of total number of signal and background events [Eq. fIlOp ]. This result 
has been used to discuss the significance of eccentricity effects and day-night variations in 
the SK, SNO, and Borexino solar neutrino experiments, as well as the sensitivity to reactor 
power variations in KamLAND. In particular, we estimate that: (1) the combination of SNO 
I+II data (306.4 + 254.2 live days) can provide indications for eccentricity effects at > 90% 
C.L., but cannot provide a determination of the day-night asymmetry accurate enough to 
constrain significantly the current LMA allowed region of parameters, even in combination 
with current SK day-night data; (2) the Borexino experiment should test eccentricity effects 
at > 3cr in three years; and (3) the KamLAND experiment should already be able to test 
time variations of the reactor power at a level > 2a. Such estimates are statistical only, and 
should be somewhat degraded by systematic uncertainties, whose evaluation is beyond the 
scope of this work. 
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APPENDIX A: ESTIMATE OF /3 ± cj/j 

In this Appendix we derive Eq. 1^. We make use of the fact that, since 



and 



N 



N 



r dt R{t) = Ns + Nb 
Jo 



Ns+Nb 



I dt 5{t-U 

•^0 ^=l 



(where 5 is a Dirac delta), we can formally write 



Ns+Nb 

dN = R{t)dt= J2 Kt-ti)dt 

i=l 



(Al) 
(A2) 

(A3) 



This "trick" allows to express the continuous function R{t) in terms of the discrete time 
series {ti}i=i^,,,^Ns+NBy to evaluate the variance of the small number R{t)dt through 
Poisson's statistics (see also 0): 



YaT{R{t)dt) = var(rfA^) = dN = R{t)dt 

Let us prove that P = Ns<Js/S: 

Ns+Nb 



E fit 
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dt 



S[l + f /(t) 
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asT 

Ns 



as- 



S 



Finally, let us prove that, for cr^/ 5 ^ 1, it is cr^ ~ y/Ns + Nb- 



var(/5) 

var 

rT 



dtR{t)f{t) 







dtf{t)R{t) 



{S + B)T + ST 

{S + B)T 
Ns + Nb 



05 

s 



dtfit) 



(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 
(AlO) 
(All) 

(A12) 
(A13) 

(AM) 

(A15) 

(A16) 
(A17) 
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In deriving Eq. ()A14|) . we have used both the distributive property var(X]i Oja^i) = 
X]iaiVar(xj) and Eq. ()A4|) . In deriving Eq. ()A16|) . we have dropped the second term on 
the right-hand-side of Eq. ()A15|) . which is not only of 0{as/S) < 1 but, in typical applica- 
tions (including all those considered in this work), is further suppressed by the small value 
of the third moment of f{t). 



APPENDIX B: RELATION WITH THE KOLMOGOROV-SMIRNOV UN- 
BINNED TEST 

The Kolmogorov-Smirnov (KS) unbinned test can be used to compare two distributions 
through the maximal difference between their cumulative distributions . In our case, the 



two distributions would be the total event rate R(t) including time variations [Eqs. (jH)) and 
(13))] and the total event rate excluding time variations, Rq = S + B. It is easy to convince 
oneself that, in general, the maximal difference d between the corresponding cumulative 
distributions is of the kind: 

where rj is a. dimensionless factor of 0(1), which depends on the functional form of f{t). 
The KS test then "measures" the difference d in units of k/\/N (for large N = Ns + Nb), 
where the number k depends on the confidence level chosen j|29|J. Therefore, the statistical 
significance of the KS test is proportional to 

dVN as Ns 

oc — — , . B2 

The comparison of the above result with Eq. (jlOj) suggests that the statistical significance 
of time variations should scale with {as/ S)/{Ns/ y/Ns + Nb), independently of the specific 
(unbinned) statistical test chosen — a fact than can be useful in prospective studies of ex- 
perimental sensitivities. Different tests may differ, however, in statistical power. In all the 
examples discussed in this work, the KS test happens to be less powerful than the one pro- 
posed through Eq. pU|) . although we cannot exclude that it might be more powerful in other 
cases. 
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